NASA Contractor Report 3389 


A Three-Dimensional Turbulent 
Compressible Subsonic Duct 
Flow Analysis for Use With 
Constructed Coordinate Systems 


R. Levy, H. McDonald, W. R. Briley, 
and J. P. Kreskovsky 


CONTRACT NAS3-21735 
APRIL 1981 




NASA Contractor Report 3389 


TECH LIBRARY KAFB. NM 

lillililllll 


□□baata 


A Three-Dimensional Turbulent 
Compressible Subsonic Duct 
Flow Analysis for Use With 
Constructed Coordinate Systems 


R. Levy, H. McDonald, W, R. Briley, 
and J. P. Kreskovsky 
Scientific Research Associates, Inc. 
Glastonbury, Connecticut 


Prepared for 

Lewis Research Center 

under Contract NAS3-21735 

(W\SA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Branch 


1981 




TABLE OF CONTENTS 





Page 

SUMMARY 1 

INTRODUCTION 2 

ANALYSIS 3 

Coordinate Systems 3 

Governing Equations .... 5 

Flow in Ducts with Superelliptic Cross Sections 9 

Potential Flow 13 

Calculations 14 

SUMMARY AND CONCLUSIONS 16 

REFERENCES 17 

FIGURES 18 

SYMBOL TABLE 28 

APPENDIX - COMPUTER CODE INPUT AND OUTPUT 29 


iii 


SUMMARY 


An approximate analysis is presented for computation of three-dimensional sub- 
sonic flow in straight and curved diffusers. The development parallels that of 
Briley and McDonald for forward-marching solution of viscous primary and secondary 
flows, but differs in the coordinate formulation used and in details of the approxi- 
mations. The present formulation is intended to facilitate the use of constructed 
coordinates in circumstances where it is difficult to maintain smooth behavior in 
higher derivatives. This analysis is applicable to nonorthogonal coordinate sys- 
tems having a curved centerline and planar transverse coordinate surfaces normal to 
the centerline. The primary flow direction is taken to coincide with the local di- 
rection of the duct centerline and is hence normal to transverse coordinate planes. 
The formulation utilizes vector components (velocity, vorticity, transport equations) 
defined in terms of local Cartesian directions aligned with the centerline tangent, 
although tliG governing equations themselves are expressed in general nonorthogonal 
coordinates. For curved centerlines, these vector quantities are redefined in new 
local Cartesian directions at each streamwise location. The use of local Cartesian 
variables and fluxes leads to governing equations which require only first deriva- 
tives of tiie coordinate transformation, and this provides for the aforementioned 
ease in using constructed coordinates. The analysis is applied to a particular 
family o^^ duct and diffuser geometries having curved centerlines and superelliptic 
cross sections, using a nonorthogonal constructed coordinate system. Present com- 
puted results are compared with available experimental measurements for different 
but related flow conditions. Qualitative agreement is observed with regard to the 
formation of secondary flows, migration of streamwise vortices, and distortion of 
the primary flow. Additional experimental data is needed to provide a basis for 
detailed quantitative evaluation of the analysis and flow predictions. 
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INTRODUCTION 


A continuing problem in the design of inlets for airbreathing propulsion sys- 
tems is the design of efficient subsonic diffusers. In the case of three-dimen- 
sional inlet diffusers, the cross-sectional shape of the ducting varies in the axial 
direction and it is frequently necessary to introduce offset bends (curved duct 
centerlines). The complicated diffuser geometry and the offset bends introduce 
strong secondary flows which have important effects on both diffuser performance 
and engine/inlet compatibility. A generalized subsonic diffuser analysis capable 
of being used as a design tool must account for several physical phenomena which 
frequently occur in practical diffusers. First, the analysis must be capable of 
treating the case when the wall boundary layers are turbulent and when the thick- 
ness is comparable to the local duct radius. Secondly, the analysis must be capable 
of treating the strong secondary flows which are commonly encountered in three- 
dimensional diffusers. 

Although solution of the three-dimensional Navier-Stokes equations would be 
definitive, more cost effective approximate methods are proving attractive. One 
such approximate method has been developed by Briley and McDonald (Ref. 1) where a 
set of governing equations are developed by introducing approximations into the 
Navier-Stokes equations and then solved by spatial marching. This paper reports 
on a study to modify the analysis of Ref. 1 for use in constructed coordinate sys- 
tems which may be nonorthogonal . By a judicious choice of dependent variables and 
of forms of the equations, some of the smoothness requirements for the coordinate 
system have been relaxed. 
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ANALYSIS 


Coordinate Systems 

A coordinate system is more than merely the grid points at which the governing 
equations are written. The coordinate system manifests itself as the geometric 
terms in the governing equations as is related to resolution of the physical pro- 
cesses and various length scales within the problem being studied. Although it is 
easy to overlook the part that the coordinate system plays in the governing equa- 
tions, the requirements for a "good" coordinate system are found by finding the 
role these geometric terms play in the fluid dynamic equations. Different forms of 
the Navier-Stokes equations can impose different requirements for a coordinate 
system. For illustrative purposes a viscous term from two forms of the two-dimen- 
sional Navier-Stokes equations will be used. 

1 2 

First consider a point P in reference Cartesian coordinates (x , x ) . The 

1 2 

generalized coordinate lines are y and y passing through P and the Cartesian 

1 2 

velocities are u and u . Using the quasi-linear form of the Navier-Stokes equa- 
tions and Cartesian velocity components one of the viscous terms in the x^ momentum 
equation would be 


j a 

Re ax^ av ^ 




where J is the Jacobian of the transformation from x^ into y^ and 3y^/3x^ is a com- 
ponent of the inverse transformation tensor. The Jacobian as well as the elements of 
the inverse transformation tensor are all related to the gradient of the position 
vector at the computational grid points relative to a fixed, reference Cartesian 
coordinate system. 
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Before continuing, let us identify this term in familiar Cartesian coordinates. 
2 2 2 2 1 

In this case J=l, dy /3x =1, y =x and u =U. By performing these simplifications 
we get: 


Re 



or 


Re dY V ^ dY / 


( 2 ) 


which is the familiar form of the dominant viscous term in the Cartesian x momentum 
equation. From the simplicity of (2) one could easily forget the complexity of the 
term when expressed in general nonorthogonal form; but the complexity is real and 
the geometric terms are present and must be properly treated. 

Returning to (1) it can be seen that the mathematical requirements on the in- 
verse transformation tensor are the same as those on the fluids variables. That is, 
each element of (1) must be continuous, have continuous derivatives and be well de- 
fined by the computational grid points being used to solve the problem. Since the 
geometric terms in (1) are directly related to first derivatives of the locations of 
the computational grid points and these are in turn differentiated, the requirement 
on both the geometry and the velocity field is that in going from one grid point to 
another the second derivatives must be continuous and well defined. In addition, 
if the discretization of the differential equations into the difference equations 
is to be accurate, the geometric derivatives must not be large compared to the fluid 
dynamic terms. 

As a second example, a corresponding viscous term is now examined in the general 

1 2 

orthogonal form of the Navier-Stokes equations. For coordinate lines y and y and 
for velocity components along the coordinate directions u^ and u^ this viscous terms 
is 

/ 

Re h I hg d Y^ I hg dY^ 



where h^ and h^ are the metrics (e.g.) 



(3b) 


To verify that this viscous term corresponds to the term examined in (1) for the 

11 2 2 

case of Cartesian coordinates the metrics are all set to unity, u^=U, x =y , x =y : 


A 



J ^ 

Re (3Y 


( 4 ) 



which corresponds to (2) , 

Examination of (3a) shows that there is a second derivative of h^ implying that 
h^ must be continuous and have continuous derivatives. Equation (3b) shows that h^ 
is itself a first derivative of the coordinate system. Thus, the requirement on the 
coordinate system for general orthogonal systems is that the coordinates must have 
continuous and well defined third derivatives. 

In constructing coordinate systems one of the constraints is the required 
smoothness. One measure of the required smoothness is the order of continuous der- 
ivatives required. The above example demonstrates that the form of the governing 
equations can directly influence the adequacy of the coordinate system. This ob- 
servation played a major role in both the selection of the forms of the governing 
equations and in the modifications made to the analysis of Ref. 1, which will now 
be discussed. 

Governing Equations 

The equations used in this analysis are based on the steady compressible 
Navier-Stokes equations. Three basic assumptions are made to permit a forward 
marching solution of the subsonic flow problem. First, the flow can be analyzed as 
being a primary flow with a transverse secondary flow. This primary flow is assumed 
to be approximately aligned with one of the coordinate directions. Second, the 
pressure field can be represented by a three-dimensional ^ priori known pressure 
field, such as comes from a three-dimensional potential flow analysis, plus a one- 
dimensional correction for viscous blockage effects and pressure loss. 

P + (5) 

The three-dimensional pressure field contains elliptic information about the flow 
(upstream influence). Third, the assumption of negligible second derivatives in the 
marching direction is also required to create an initial boundary value problem which 
can be solved by forward marching in one of the coordinate directions. 

With these assumptions the Navier-Stokes equations are written using Cartesian 
velocities, u^, and Cartesian coordinates, x^. Although subsequent forms of these 
equations are to be evaluated in general nonorthogonal coordinates, y , the use of 
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Cartesian velocities and Cartesian directions for the vector equations will result 

in the favorable smoothness properties discussed above. The Cartesian coordinate 

system x^, is rotated at each computational step to keep the Cartesian direction 

aligned with the tangent to the duct centerline ensuring that even when the duct 

3 

contains large bends, the x coordinate direction is always roughly aligned with the 
primary flow direction. The x^ coordinate direction is aligned with the outward 
facing normal. Special treatment of cases where the centerline curvature is zero 
is readily provided. The vector quantities, velocity and vorticity, are also 
transformed to stay consistent with the Cartesian directions. 


The velocity vector, V, is decomposed into a primary flow component, u^, normal 
to the cross plane and a secondary flow component. The secondary flow is further 
decomposed into parts associated with the irrotational cross flow, V , , and the ro- 
tational cross flow, V,: 




V = u^ + Vs 


V = 0, + + V 




( 6 ) 

(7) 


where 


V3X( v^) = 0 


Vs-(pV^) = 0 


( 8 ) 

(9) 


and is the surface gradient operator in the cross plane: 

„ d T d T dr^ d V dY^ d T 

^ i ■ -.wi -,v,n ' -,^^2 j 


ax' 


ax‘ 


ax' ar^ 


ax*' av' 


( 10 ) 


for n=l,2. 

3 3 

Since the x direction is the primary flow direction, the x momentum equation 

in Cartesian coordinates is 


du- 


) u . 


j axJ 


ap,_ 


3-D 


a X' 


ap 


I- D 


ax' 


-t- F, 


axJ 


Me 


- s 




3 , 

j ax^ 


auk 

3 


^ Sjjd-SK^) 


( 11 ) 
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where is the body force in the x direction and u is the effective viscosity. 

J 63 

Using the nonorthogonal computational coordinates, y^, the x momentum equation: 




)U- 


j axJ ayi 


ay' aP3- D 
ax^ av ' 


av* ap 


l-D 


ax^ av' 


+ F, 


ay' a 


axJ ay 






j/\axJ ay^ ax^ ay" 


( 12 ) 


-2. » ( I _ R , ) auk 

3 3j k ijyk ^yn 


is solved along with the global mass flux conservation relation 


the energy equation: 


II ‘ ° 


= const. 


(13) 


(14) 


and the equation of state: 


P = Pg.o + P,_o = ^RT 


(15) 


to determine U 3 , p, and )■ 

The continuity equation is written as: 

+ 73 -(/ 3 V^) - 0 

Introducing a potential function cj): 
this Poisson equation is solved in the form: 


du^p 


(16) 


(17) 


ayJ a 


ax' ayJ V ax' / ayk 


d(f> 

+ nl 

^ay'' 


ayk 

+ P\ 

^ ax* 

) 7^ 


ay' 

ay2 

d^4> 

+ 2p 

ax' 

ax' 

ay'ay ^ 


ay" dpuj, 
ax^ ay" 


( 18 ) 
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( 19 ) 


for k=1^2 where <(> is defined by: 
and satisfies (8) . 

The remaining unknown is the rotational cross flow velocity component, V . 

12 ^ 

This is found by cross differentiating and combining x and x momentum equations 

3 

to form the vorticity transport equation in the x direction. One form of this 
equation is: 




- n 


d u 


dp rij 


J dX^ 


J axJ 


Up d^Sl^ 

T 


a(F,//D) a(F2/^) 


ax' 


ax' 


p 


a/Xe a^e 


ax' 


dX‘ 


( d\i\ augN auj auj 2 




ax' 




+ 

P 


axJ 


,ax' 


au I 

a^ 


2 auk 


a / duz 

Jx' \ a^ 


3 axk 


( 20 ) 




du\ 


du j 


duk 


p 

ax^axj 

V axj 

I 0|i 

ax' 3 'J 

dX^ 

\ 


/ au2 

, <^^j 2 , 

auk 

p 

ax' axj 

\ axj 

ax2 3 -2j 

dX^ 

1 

/ dp 

dPj-D 

dp dP3-D \ 



\ dX^ 

ax' 

ax' ax^ ) 



where u^ and u^ are components of 


u, i + Ugj = Vj = 


( 21 ) 
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and the vorticity is defined as: 


^2 

To solve this equation for V, 
introduced satisfying: 


dU2 

auj 


dX^ 

ax® 

(22) 

dU3 

dU| 


1 

X 

II 

ax^ 

(23) 

du I 

(3u 2 



ax' 

(24) 
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a stream function, i/j in the x direction, is 


^ 5 ^ 4 ' = 


(25) 


Using (24) the stream function and the nominally streamwise component of vorticity, 
are related by 


n^^VsXVs 

= Vg X ( ) 

(26) 

= VgX 

^3 = Vs X Vs X 


Equations (20) and (26) are solved simultaneously for ijj and Since the remaining 

parameters are either known functions or expressible in terms of known functions, 
and ij). 

Fl ow in Ducts w ith Supere ll iptic Cross Sections 

The analysis of the previous section was applied to the flow in a general class 
of subsonic diffuser ducts. The duct centerline is defined by polynomials and the 
cross sections are described by superellipses in a plane locally normal to the cen- 
terline, The duct boundaries are coordinate surfaces, greatly simplifying the 
boundary conditions. This coordinate system is described in detail in Ref. 2. 

Note that the grid distribution functions of Ref. 2 have been replaced. For radial 
distributions to resolve the boundary layer: 
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R 


tonh(Dr) 

tanh(D) 


(27) 


where R is physical radius, r is the radial computational coordinate parameter and 
D, the damping function which determines the severity of the distribution, is con- 
trolled through the parameter a by: 

“ = «« 

A circumferential distribution function is also used to cluster grid points 
into regions of tight wall curvatures. The equally spaced computational coordinate 
t is related to the physical angle 0 in a quadrant by the relation: 


® = A sinh 


/p(t-c) \ 


+ B 


t, =0, tg= I 

@ 1=0 @2 = tt /2 

®j.= tan"' (S^) 

S = shape (ratio of major to minor axis) 
^ - superelliptic exponent 

®| - @r 

R= 7 T T- , D = In (2 3KE + 23KS) 

®2-®c 

ln(2*E +2*S) 


(29) 


where : 


DZM 



sinh(DZM-t|) - R» sinh(DZM-t2) -i 
cosh(DZM- 1 |) - R • cosh(DZM-t2) -I 


(30) 


@2“ ® I 

sinh (DZM -Ctg -C)) (31) 

sinh ( DZM-(t| -O) 
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©2 - ®l 



sinh (DZM(t 2 -C)) “ sinh (DZM(t, - O) (32) 

This distribution function, (29), clusters grid points toward 0=©^- Several examples 
of the effects of (27) and (29) are presented in Fig. 1 for a series of cross sec- 
tional shapes. 

The governing equations used in this application are: 

3 

MOMENTUM in the X direction 


av' auj 

ayi 


av‘ aP3-D _ ^ aP|_D 
ax^ av ' ax^ ar^ 


ay' 

a 

r f( 


au3 

ay" 

auj \ 


axl 

ay' 


1 - Si3)( r 

v\axi 

— ^ + 
ay" 

ax^ 

ay" j 

(33) 


-f 




av" ^ 
^x’' ay" 


from CONTINUITY 

ay^ a / a^\ a</> ^ /ay*'\^ d^ 4 > 

ax' ayJ v ax' / ay*' ^ ax' / Jy*^ 

ay' ay 2 d^4> 
ax' ax' ay'ay^ 


(34) 

ay*^ dpu^ 

ax^ ay" 


VORTICITY TRANSPORT in the X direction 

a U 


a ^3 

axi 


- n 


J axJ 


STREAM FUNCTION 


H-e a^fls 

P axj 2 


ay^ a 
ax' ay) 


I av*^ \ ai// I / ay*' a^V' 
7 ax' / "ay^ ^ \ ax' / ay**^ 


2 ay' ay^ a^</^ 
ax' ax' ay'ay^ 


(35) 


(36) 
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The vorticity transport equation was used in simplified form for convenience during 
the development phase of this study. Use of equation (20) for vorticity transport 
is anticipated. 

At each marching step the computational domain has two symmetry boundaries, a 

boundary which is the duct wall and a boundary which is a small tube surrounding 

3 

the centerline singularity (see Fi. 2). The boundary conditions for the X momentum 
equation are ^^u/9n=0 on the two symmetry surfaces where n is the local surface 
normal, u = 0 on the wall and 

du _ _ J_ dP 

^ ' P 

at the centerline. 

The scalar potential equation for 4> has the boundary condition (j)=0 on all 
boundaries . 

The vorticity and stream function equations are solved sequentially and coupled 
iteratively. The stream function boundary condition is applied at all bound- 

aries. The vorticity boundary condition ^^=0 is applied along the symmetry and 
centerline boundaries. On the wall the vorticity boundary condition is derived fol- 
lowing Briley & McDonald (Ref. 1), by writing the stream function equation (36) on 
the wall using central difference formulas. The additional constraint needed to 
eliminate the added variable is the no-slip condition for the cross flow. These 
equations are 


ayj 

d 

(± 

ay2 \ 



ax' 

ayj 


axi / 

ay^ P ' 

^ ax' ) 

2 

ay' 

iL" 


iz. 3 


P 

ax' 

ax' 

ay'ay^ 



(37) 


and setting the slip velocity from the solution of the stream function-vorticity 
equations equal and opposite to the slip velocity from the solution of the scalar 
potential equation yields: 



P 


av' 


dx' 

ax2 


ayf 

ax' 


ay' 

ax^ 

. ^ 
’ ax' 


ax 


2 

2 


(38) 
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where 


^ V ' vi 


(39) 


9 ' ax' av^ 9 ' ax^ av^ 

on the wall. 

The mixing length model used in this analysis employs the eddy-viscosity formu- 
lation for the Reynolds stresses, i.e. , 

My f3v i 

axi 

Hence, this formulation still suffers from the physical shortcoming that there 
is zero Reynolds stress wherever the velocity gradient is zero. In addition, the 
eddy viscosity formulation is isotropic which may be incorrect in many three-dimen- 
sional and swirling flows. However, for practical calculations of turbulent inter- 
nal flows there are as yet no other transport models which are either suitable or as 
relatively well developed. 

The mixing length turbulence model employed in this analysis is based on a 
mixing length distribution. The mathematical form of the expression for the tur- 
bulent viscosity follows from Ref. 3: 


= = 1/2 

— ; pi.2(2e:e) 

where e is the mean flow rate of strain tensor 

e= 1/2 [(Vv) +(Vv)''’] 

The mixing length Z is determined from the Buleev formula (Ref. 4) 


I 


2k 


1 ^ dS 


(40a) 


(40b) 


(41) 


where K is the von Karman constant and S and 0 are defined in Fig. 3. 


Potential Flow 

The present application of the analysis presented in the previous section uses 
a three-dimensional potential flow to supply the _a priori known pressure field used 
in the viscous flow analysis. The potential flow is computed in the computational 
coordinate system by solving the equation 

= 0 (42) 


d^4> d^4> d^4> 

dx'^ ^ ax2^ ^ dx^^ 


0 


or expressed in the computational coordinate system: 


(43) 
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( 44 ) 


where 4) is the velocity potential, are the Cartesian coordinates and are the 

computational coordinates. This equation is solved using an iterative scalar ADI 

procedure. The boundary conditions at the duct wall is that there is no flow 

d d> 

through the wall, i.e. = 0 • Since the potential, (}), is dimensionless 

d n 

the inflow and outflow values are set to 0 and 1 respectively. The duct centerline 
would represent a coordinate singularity so the coordinate system used extends from 
the duct walls down to a small tube of finite radius around the centerline. The 
boundary condition applied at the centerline is derived from a zero velocity normal 
to the symmetry plane and permitting a velocity component tangent to the symmetry 
plane. Since the solution technique is iterative, the velocity boundary condition 
is lagged and updated every few iterations. The area in the tube around the center- 
line is small and this explicit, rather than implicit, handling of the boundary 
condition appears to have little if any adverse influence on the observed conver- 
gence properties of the solution procedure for solving Eq. (44). 

Once the converged value of the potential, cj) , is found, the non-dimensional 
velocities at each grid point are computed. These velocities are scaled to be 
compatible with the specified reference conditions- From the velocities the cor- 
responding pressure field is computed and values of pressure coefficient are com- 
puted and saved on permanent storage for later use in the viscous spatial marching 
calculations . 

Calculations 

The above analysis was applied to a series of test problems. The first case 
presented here is the flow in the entrance region of a straight circular pipe. The 
calculations were compared to the experimental data of Reshotko (Ref. 5). The ex- 
cellent agreement with the data is shown in Fig. 4 . 
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The second case was of turbulent flow through a circular pipe undergoing an 
S-shaped bend. The bend was offset 1 diameter in a streamwise distance of 5 dia- 
meters (Fig- 5). The entrance Mach number was 0.2 and the entrance Reynolds number 
was 10^. Although no test data was available to quantitatively evaluate the re- 
sults, the test case was apparently successful, A qualitative evaluation has been 
made by comparing salient features of the computed flow with the flow measured in a 
model of the Boeing 727 Center Duct Inlet (Ref. 6). It must be emphasized that the 
computed case was for different geometry, Reynolds number and Mach number than the 
tested case. However, not only are the classical secondary flow patterns clearly 
discernible in both cases, but the development of the secondary flow patterns ap- 
pears to occur in the same manner. The locations marked A, B, C and D refer to 
cross sections where comparisons will be presented. Sections A are near the center 
of the first bend, sections B are just after the inflection, sections C are just 
after the center of the second bend, and sections D are at the ends of the ducts. 
Figure 6 presents the contours of total pressure at the bottom of the 727 duct for 
stations A, B and C. The growing region of low total pressure is clearly seen. 

Also shown are contours of streamwise velocity at the bottom of the computed S-shaped 
duct. In the case presented these streamwise velocity contours have nearly the 
same shape as total pressure contours. Here again the growing region of low total 
pressure is clearly seen at the bottom of the duct. A comparison of sections D in 
Figure 7 shows the classical total pressure contour shapes for S-shaped bends. 

Figure 8 also shows these classical exit plane total pressure contours from three 
other S-shaped bends (Ref, 7). 
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SUMMARY AND CONCLUSIONS 


An analysis is presented for computing three-dimensional subsonic flow in dif- 
fusers. Following Briley and McDonald, approximations are introduced into the 
Navier-Stokes equations to produce a system of equations which can be solved by 
forward marching. This paper reports on a study to adapt the above analysis for 
use with nonorthogonal constructed coordinate systems. 

At each computational plane a local Cartesian coordinate system is defined at 
the duct centerline with one coordinate direction tangent to the centerline and 
another coordinate direction aligned with the outward facing normal. Special treat- 
ment of cases where the centerline curvature is zero is readily provided. Velocity 
components are defined to be aligned with these Cartesian directions and those 
governing equations that are vector equations are written in these local Cartesian 
directions. At each marching station the vector quantities computed at the previous 
station are transformed into these new vector directions. 

The resulting partial differential equations are transformed into the construc- 
ted coordinate system. These equations contain only first derivatives of the ele- 
ments of the transformation from Cartesian to constructed coordinates as contrasted 
with the familiar general orthogonal equations which contain second derivatives of 
the metric scale factors (h’s). This reduced sensitivity of the governing equations 
to the details of the coordinate system permits more flexibility in the coordinate 
construction process. A family of duct shapes with superelliptic cross sections is 
currently available within the code representing a large class of diffuser shapes. 
These coordinates are nonorthogonal. The analysis developed in this study is ap- 
plicable to more general cross-sectional shapes than are currently available within 
the code. 

The above analysis was applied to subsonic viscous flow in several ducts. 
Qualitative agreement with experimental measurements was encouraging since the 
formation and migration of the vortex seen in the measurements was also present in 
the computed results. Additional experimental data is needed to provide a basis for 
detailed quantitative evaluation of the analysis. 

The current version of the computer code solves the streamwise vorticity and 
stream function equations by coupling them through a sequential iterative procedure. 
A preferred method would be to solve these equations as an implicitly coupled block 
system. This is expected to improve convergence properties of these equations and 
permit larger streamwise steps to be taken. 
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SYMBOL TABLE 


y 

i, j, k 

'ij 

V 


u, 

1 


p 

p 

Re 

y 

y^ 


1 


r 


t 


Cartesian coordinate 

General curvelinear coordinate 

Cartesian unit vectors 

Kroniker delta 

Velocity vector 

Velocity component 

Pressure 

Density 

Reynolds number 
Viscosity 

Effective turbulent viscosity 

Vorticity component 

Scalar potential function 

Stream function 

Radial computational 
coordinate parameter 

Azimuthal computational 
coordinate parameter 
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APPENDIX 


COMPUTER CODE INPUT AND OUTPUT 


The PEPSIC code input is of the form: 
MODE, TITLE 
Namelist Data 


MODE, TITLE 
Namelist Data 


MODE, TITLE 
Namelist Data 


STOP 

The mode parameter selects which of the four functions available within the 
code are to be used. 

MODE - 4 geometric preprocessing 

3 geometry written for potential flow on logical unit #10 
2 potential flow calculation 
1 parabolized Navier-Stokes 
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\<fhen operating in MODE = 4, the code will solve a series of simultaneous 
equations to determine polynomials Y = f (X) from tabulations of Y, Y* and/or 
Y” versus X. The input is specified through NAMELIST: 

NEQ - number of equations 
XIN - input table of X values 

YIN - input of Y, Y^ or Y" according to corresponding value of ID 

ID - 0 YIN is Y 

1 YIN is Y’ 

2 YIN is Y” 

NX - number of locations in the interval XI, X2 at which the resulting polynomial 
is to be evaluated. 

XI, X2 - interval in which the resulting polynomial is to be evaluated 

The input for MODE = 3, 2 and 1 is provided through three Namelist tables. 
Although the code is still in the development stage and elegant input has not yet 
been devised, the Namelist method provides a serviceable means of setting up the 
cases to be run. The three Namelists are, in order of being read: GEOM, FLUIDS, 

and SWITCH. The numbers in parenthesis following the explanation of the input 
variable are the recommended values. 


$GEOM 

NRAD - Number of grid points in the radial direction (21) 

MTHETA - Number of grid points in the circumferential direction (21) 

RTUBE - Radius of the tube put around the centerline to remove the 

geometric singularity (0.01) 

P(l, J) - Coefficients of polynomial in x for shape factor - ratio of major to 
minor axis where: 

S.F. = P(l,l) + P(l,2) X + P(l,3) + P(l,4) x^ . . . 

+ Pd.ll) 


P(2, J) - 
P(A, J) - 

PA(3,1,J)- 


Coefficients of polynomial in x for radius of major axis 
Coefficients of polynomial in x for superelliptic exponent 
between 2 and 10 

Coefficients of polynomial in x for transverse Cartesian component 
of centerline curve. Set to zero for straight duct. 
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NOTE: The P and PA arrays determine the location of the grid points. Since the 

position vector of the grid points is differentiated twice within the code, the 
piecewise polynomials describing these locations must be continuous to two 
derivatives. Small values of higher derivatives is also recommended. 

YZERO - Reference length in ft. (1.0) 

NS - Number of the last marching station to be computed 

NFIRST ~ Number of the first marching station whose geometry is to be 

specified (usually 1) 

X (NFIRST) - Nondimens ionali zed location of first marching step whose geometry 
is to be specified 
DX - Initial step size 

AP - Ratio of consecutive step sizes in the marching direction 

x(i) = x(i-l) + AP* [x(i-l)-x(i-2) ] 

IPA - Circumferential distribution function, 2 is equally spaced, 

3 is automatic distribution, (3) 

VIS - a in Eq. (28) (.3) 

IPB “ Radial distribution function 2 is equally spaced, 3 is packed toward 

the wall (3) 

IRSTIN - The station number to be read in for restart 
If IRSTIN = 0, not a restart case 
IRSTOT - The interval for saving restart information 
If IRSTOT = 0, no restart information saved 
JRSTIN - Logical file name from which restart information is to be read 

JRSTOT - Logical file name onto which restart information is to be written 

NFILE - The sequence number in JRSTIN of the desired restart information 

NSAVED - The number of restart stations saved on JRSTOT; must be initialized in 

inputs to the number of stations already written (and to be preserved) 
on JRSTOT. Nominally initialized NSAVED = 0. 

NOTE: By setting JRSTOT = JRSTIN and NFILE = NSAVED, one file can be used for both 

read and write information, without destroying the already saved information. 
$FLUIDS 

3 

RZERO - Reference density in Ibm/ft 

UZERO - Reference velocity in ft/sec 

CMACH “ Reference Mach Number 
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NOTE: Set either CMACH or UZERO and the other will be computed within the code. 

REY “ Reynolds number based on UZERO, YZERO and RZERO, Reference viscosity is 
determined by VISCOS = UZERO*YZERO+RZERO/REY 
KTURB - 0 Laminar 

- 1 Turbulent 

YSL0T(2)“ Initial nondimensionalized boundary layer thickness 
$ SWITCH 

NEIPF - Number of circumferential grid points used in potential flow 
NE2PF - Number of radial grid points used in potential flow 
IGDMP “ 2 Dumps strearawise momentum equation coefficients loaded into 
matrix inverter 
0 No dump 

IPLOT - 1 Write plot file (logical unit //8) 

0 No plot file 

IC0EF(1,N) - Auxiliary array which controls various options within the code. 

In general a value of 0 means the option is off, a value of 1 means on. 
N 

1 COEFG convective terms printout 

2 COEFG first viscous terms printout 

3 COEFG second viscous terms printout 

4 COEFG inviscid pressure gradient printout 

5 COEFG viscous pressure gradient printout 

6 PFIELD, CPI=0 

7 PFLOW inter dump -1, full dump 1 or 2 

8 El dump from CROSEC 

9 output PHI, PSI velocities 

10 output PHI, PSI Cartesian velocities 

11 output Cartesian velocities 

12 ADI dump PSI 

13 ADI dump vorticity 

14 ADI dump PHI 
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15 input NAMELIST dump 

16 IPRFSE dump 

17 DYFA mixed der. diagnostic message 

18 level 2 output 

19 PFIELD drum reading information 

20 output for ITER=1,2 

IC0EF(2,N) N 

1 wall vorticity for no slip (COEFVS) 

2 No wall damping on NLEN 

3 COEFVS iteration print 

4 Buleev mixing length formula 

5 number of pressure iterations in ADICUP 

6 auxiliary variable printout 

7 auxiliary variable printout 

8 auxiliary variable printout 

9 auxiliary variable printout 

10 auxiliary variable printout 

11 auxiliary variable printout 

12 auxiliary variable printout 

13 auxiliary variable printout 

14 PFLOW intermediate solutions 

15 FRAME printout 

16 ADI convergence information 

17 DPDXV iteration printout 

18 printout 3 level velocity and vorticity 


NOTE: 

It is advisable when running the potential flow calculation for' the elliptic 
pressure field (M0DE=3,2) to add constant area extensions both upstream and doxvTi- 
stream of the duct to be evaluated in order to minimize the effects of the upstream 
and dovmstream boundary conditions. 
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The output for M0DE=4 presents the coefficients of the polynomial and an 
evaluation of the polynomial at NX points in the range from XI to X2. Displayed 
are the values of the polynomial along with the first and second derivative. 


NEQ= 6 
Xl= 0.0 
X2= 10.0 

Hx= n 

XIH= 3X0.0, 3X10.0, 2X0.0 
YIM= 3X0.0, 2.0, <iX0.0 
ID= 0, 1, 2, 0, 1, 2, 2X0 
4EHD 
ANSWER 
ANSWER 
ANSWER 
ANSWER 
ANSWER 
ANSWER 


1 
1 

2 

3 

4 

5 

6 
7 
S 
9 

10 

II 


0.12000000E-03 
-0.30000000E-02 
Q.20000000E-01 
0 .00000000 
0.00000000 
0.00000000 


K XXX 
X XXX 
X XXX 
X XXX 
X XXX 
X XXX 


0.00000000 
O.IOOOOQOOE 01 
O.2OOOQO0OE 01 
0.30000000E 01 
0.90000000E 01 
0.50000000E 01 
Q.60000000E 01 
0.70000000E 01 
O.aOOOQOQOE 01 
0.90000000E 01 
O.IOOOOOOOE 02 


0.00000000 
0.17120000E-01 
O.X1589000E 00 
0.326I6000E 00 
0.63988000E 00 
O.IOOOOOOOE 01 
0.13651200E 01 
0.16738900E 01 
0.188916Q0E 01 
0.19828800E 01 
0.20000000E 01 


YP 

0.00000000 
Q.98600000E-01 
Q.15360000E 00 
0.26960000E 00 
0.3956Q000E 00 
0.37500000E 00 
0.39560000E 00 
0.26960000E 00 
0.15360000E 00 
0.98600000E-01 
“0.888178926-15 


YPP 

0.00000000 
0.86900000E-01 
0.I1520000E 00 
0.10080000E 00 
0.57600000E-01 
0.27755576E-16 
-0.57600000E-01 
-0.10Q80000E 00 
-0.11520000E 00 
-0.86900000E-01 
“0.66613381E-15 


The output for MODE=3 provides summary tables of geometric information includ- 


ing 

the absolute value 

of the radius 

of curvature of 

the centerline. 


JX 

X 

SHAPE 

MAJOR AXIS 

ROTATION 

EXPONENT 

9 

0.00000000 

O.IOOOOOOOE 01 

O.IOOOOOOOE 01 

0.00000000 

0.20000000E 01 

5 

0.99000000E 00 

O.IOOOOOOOE 01 

O.IOOOOOOOE 01 

0.00000000 

0.20000000E 01 

6 

0.1980Q000E 01 

O.IOOOOOOOE 01 

O.IOOOOOOOE 01 

0.00000000 

0.20000000E 01 

7 

0.29700000E 01 

O.IOOOOOOOE 01 

O.IOOOOOOOE 01 

0 .00000000 

0.20000000E 01 

8 

0.39600000E 01 

O.IOOOOOOOE 01 

O.IOOOOOOOE 01 

0.00000000 

0.20000000E 01 

9 

0.99500000E 01 

O.IOOOOOOOE 01 

Q.lQOOQQOOE 01 

0.00000000 

0.20000000E 01 

10 

0.59900000E 01 

O.IOOOOOOOE 01 

O.IOOOOOOOE 01 

0.00000000 

0.20000000E 01 

11 

0.69300000E 01 

O.IOOOOOOOE 01 

O.IOOOOOOOE 01 

0.00000000 

0.20000000E 01 

12 

0.79200000E 01 

O.IOOOOOOOE 01 

O.IOOOOOOOE 01 

0.00000000 

0.20000000E 01 

13 

0.89100000E 01 

O.IOOOOOOOE 01 

O.IOOOOOOOE 01 

0.00000000 

Q.20000000E 01 

19 

0.99000000E 01 

O.IOOOOOOOE 01 

O.IOOOOOOOE 01 

0.00000000 

0.200000QQE 01 




CENTERLINE LOCATION 

(X1,X2,X3) 



JX 

X 

XI 

X2 

X3 

RAD OF CURVE 

9 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

5 

0.99000000E 00 

0.16638311E-01 

0.00000000 

0.99000000E 00 

0.11688730E 02 

6 

O.I9800000E 01 

0.11279103E 00 

0.00000000 

0.19800000E 01 

0.89889780E 01 

7 

0.29700000E 01 

0.31826750E 00 

0.00000000 

0.29700000E 01 

0.10856523E 02 

8 

0.39600000E 01 

0.62110269E 00 

0.00000000 

0.39600000E 01 

0.19796253E 02 

9 

0.99500000E 01 

0.98125125E 00 

0.00000000 

0.99500000E 01 

0.90607325E 03 

10 

0.59900000E 01 

0.13992822E 01 

0.00000000 

0.59900000E 01 

0.21S37661E 02 

11 

0.69500000E 01 

0.16550729E 01 

0.00000000 

0.69300000E 01 

0.112906Q9E 02 

12 

0.79200000E 01 

0.18715030E 01 

0.00000000 

0.79200000E 01 

0.90Q87385E 01 

13 

0.89100000E 01 

0.19781995E 01 

0.00000000 

0.89I00000E 01 

fl.ll025311E 02 

19 

0.99000000E 01 

0.19999803E 01 

0.00000000 

0.99000000E 01 

0.85S92988E 02 
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upon successful completion of writing the geometry on file #10 the following 
message is written: 

GEOMETRY WRITTEN FROM STATION 5 AT X=0.9900E 00 TO STATION 14 AT X=0.9900E 01 FOR ITCNT= 5 


The output for MODE=2 consists of intermediate printout information about the 
iteration followed by a tabulation of the information written on the pressure field 
file logical unit 13. Listed at each transverse computational plane are the stream- 
wise plane number, the parametric location of the centerline, the radius from the 
centerline of each row of grid points expressed as a fraction of the local duct 
radius, the azimuthal location in radians of each radial line of grid points, and a 
two-dimensional display of pressure coefficient, Cp, at each node used in the poten- 
tial flow calculation. 


ADI ITER= 1 PMAX= O.IOOOQE 01 DELP= 0.115S<iE-03 RH5MAX= 0.20776E-O2 DUI11= 0 . '•'iS 15E-0 1 DUM2= 0.11151E 01 

ADI ITER= 2 PMAX= O.IOOOOE 01 DELP= 0.3A012E-0'* RHSMAX= 0.55896E-03 DUM1= 0 . 16<*8E- 01 DUt12= 0.10323E 01 

ADI »» ITERS 3 PMAX= O.IOOOOE 01 DELP= 0.10552E-0‘* RHSnAX= 0.20576E-03 DUni= 0 . 1 E-0 1 DUN2= 0.10257E 01 

ADI AK ITER= <t PMAX= O.IOOOOE 01 DELP= C.3Za67E-05 RH5MAX= 0.6445^E-04 0UM1= O.-iOSOZE-Ol DUM2= 0.10Z99E 01 

ADI «« ITER= 5 PMAX= O.IOOOOE 01 DELP= 0.10287E-05 RHSMAX= 0.20220E-OA DUMls O.A0751E-01 DUM2= 0.10175E 01 

ADI ITER= 6 PMAX= O.IOOOOE 01 DELP= 0.10003E-03 RH5MAX= 0.2Q8A2E-02 DUM1= O.Sa'VAOE-Ol DUM2= 0.95987E 00 

ADI ITER= 7 PMAX= O.IOOOOE 01 DELP= 0.30161E-0A RHSMAX= 0.660A5E-03 DUMls 0.36932E-01 DUM2= 0.91333E 00 

ADI wx ITER= 8 PMAX= O.IOOOOE 01 DELP= 0.93756E-05 RHSMAX= 0.20635E-03 DUM1= 0.36582E-01 DUM2= 0.90871E 00 

ADI KX ITER= 9 PMAX= O.IOOOOE 01 DELP= 0.29228E-05 RH5MAX= 0.6'i658E-0A DUM1= 0.36338E-01 DUM2= 0.90A06E 00 

NONLINEAR INNER CONVERGENCE IN 10 ITERATIONS 
ADI IT0UT= 6 DELPQ= 0.29375E-Q3 DUM1= 0.29373E 00 

NONLINEAR ELLIPTIC EQUATION CONVERGES IN 6 ITERATIONS 


SOLUTION AT LZ= 1 1-0.29700E 01 

0.00000 O.lllll 0.22222 ' 0.33333 O.-iAAA-V 0.55556 0.66667 0.77778 0.88889 1.00000 

0.83267D-16 0.38638 0.70867 1.0186 1.3661 1.7755 2.1230 2. ‘♦329 2.7552 3.1A16 

10 0.83228E-03 0.80667E-03 0.7298AE-03 0.6218AE-03 O.A7708E-03 0.30190E-03 0.15920E-03 0 . 58856E-OA-0 . 67039E-05-0 . 28558E-0<i 
9 0,82123E-03 0.79608E-03 0.72069E-03 0 . 6 l^tSA E~03 0.9722AE-03 0.50001E-03 0.15962E-03 0 . 60859E-0‘i-0 . 36895E-05-0 . 25206E-OA 
8 0.78807E-03 0.76A32E-03 0.69306E-03 0.59265E-03 O.A5772E-03 0.29A31E-03 0.16088E-Q3 Q-SOSOVE-O^i 0 . 53539E-05-0 . 15152E-0A 

7 0.73351E-03 0.71197E-03 0.6't736E-03 0.55610E-03 O.A3313E-03 0.28398E-03 0.16183E-03 0.75570E-OA 0 . 1 905A E-O't 0.21A60E-06 

6 0.65a60E~03 0.63997E-03 0.58‘iOaE-03 0.50996E-03 0.39799E-03 0.26789E-03 0.16087E-03 0.85077E-09 0.3536‘iE-04 0.1S793E-04 

5 0.56A39E-03 0.5^921E-03 0.50369E-03 0 .A3916E-03 0.35158E-03 O.a'iASAE-OS 0.15595E-03 0.92993E-0A 0.51697E-0‘i 0.37932E-04 

4 0.45153E-03 0.44019E-03 0.40616E-03 0.35820E-03 0.29285E-03 0.21215E-03 0.14460E-03 0.96441E-04 0.65055E-04 0.54593E-04 

3 0.31981E-03 0.31250E-03 0.29057E-03 0.26082E-03 0.22055E-03 0.16892E-03 0.12418E-03 0.92483E-04 0.72659E-04 0.66051E-04 

2 0.16733E-03 0.16441E-03 0.15566E-03 0.14635E-D3 0.13504E-03 0.11597E-03 0.96157E-04 0.83313E-04 0.77106E-04 0.75038E-04 

1 0.00000 0.53068E-05 0.21227E-04 0.49795E-Q4 0.86259E-04 0.11860E-03 0.14521E-03 0.16972E-Q3 0.18833E-03 0.19453E~03 


SOLUTION AT 

II 

ro 

Z-0.19800E 

01 







0. 00000 

0.11111 

0.22222 

0.33333 

0.44444 

0.55556 

0.66667 

0.77778 

0.88889 

1 . 00000 

0. 832670-16 

0.38638 

0.70867 

1.0186 

1.3661 

1 .7755 

2.1230 

2.4329 

2.7552 

3.1416 


10 0.61894E-03 0.70333E-03 0.95646E-03 0.13282E-02 0.18386E-02 0.24482E-02 0.29393E-02 0.32925E-02 0.35335E-02 0.36139E-02 

9 0.64007E-03 0.72310E-03 0.97214E-03 0.13580E-02 0.18403E-02 0.24403E-02 0.29258E-02 0.52715E-02 0.35089E-02 0.35880E-02 

8 0.70346E-03 0.78241E-03 0.10192E-02 0.13672E-02 0.18452E-02 0.24165E-02 0.28773E-02 0.32088E-02 0.34351E-02 0.35105E-02 

7 0.80568E-03 0.87800E-03 0.10950E-02 0.14140E-02 0.18525E-02 0.23769E-02 0.28007E-02 0.31057E-02 0.35138E-02 0.33831E-02 

6 0.94109E-03 0.10046E-02 0.11949E-02 0.1475ZE-02 0.18607E-02 0.23221E-02 0.26956E-02 0.29646E-02 0.31482E-02 0.32093E-02 

5 0.11022E-0Z 0.11549E-02 D.13132E-02 0.15A6AE-02 0.18676E-02 0.22523E-02 0.25644E-02 0.27892E-02 0.Z94Z7E-O2 0.29939E-02 

4 0.12795E-02 0.13202E-02 0.14421E-02 0-16222E-02 0.18705E-02 0.21681E-02 0.24097E-02 0.Z5838E-02 0.27030E-02 0.27427E-02 

3 0.14617E-02 0.14893E-02 0.15720E-02 0.16955E-02 0.18664E-02 0.20701E-02 0.22348E-02 0.23539E-02 0.24362E-02 0.24636E-02 

2 0.16341E-02 0.16486E-02 0.16918E-02 0.17586E-02 0.18527E-02 0.19622E-02 0.20481E-02 0.21112E-02 0.21562E-02 0.21712E-02 

1 0.17884E-02 0.17929E-02 0.18052E-02 0.18284E-D2 0.18652E-OZ 0.19096E-02 0.19376E-02 0.19561E-02 0.19718E-02 0.19770E-02 


SOLUTION AT 
0.00000 

LZ= 3 

O.lllll 

3-0 . 99000E 00 
0.22222 0.33333 

0.44444 

0.55556 

0.66667 

0.77778 

0.88889 

1 .00000 

0.83267D-16 

0.38638 

0.70867 1.0186 

1 . 3661 

1.7755 

2.1230 

2.4329 

2.7552 

3.1416 


10-0.35456E-02-0.30759E-02-0.16674E-02 0.35862E-03 0.30963E-02 0.63416E-02 0.89146E-02 0.10742E-Q1 0.11978E-01 0.12390E-01 

9-0.35965E-02-0.29358E-02-0.15546E-02 0.43366E~03 0.31218E-02 0.63094E-02 0.88387E-02 0.10636E-01 0.11852E-01 0.12257E-01 

8-0.29492E-02-0 .25156E-02~0. 12155E-02 0.65960E-03 0.31992E-02 0.62135E-02 0.86117E-02 0.10318E-01 O.U472E'01 0.11856E-01 

7-0.22555E-02-0. 18635E-02-0.68812E-03 0.10118E-02 0.33197E-02 0.60631E-02 0.82542E-02 0.98154E-02 0.10872E-01 0.112Z4E-01 

6-0.136A5E-02-0.10254E-02-0.87377E-05 0.14662E-02 0.34745E-02 0.58668E-02 0.77859E-02 0.91561E-02 0.10084E-01 0.10392E-01 

5-0. 32330E-03-0.45531E-04 0.78728E-03 0.19994E-02 0.3655ZE-02 0.56323E-02 0.72261E-02 0.83657E-OZ Q.91392E-02 Q.93966E-02 

4 0.8Z2O7E-O3 0.10329E-02 0.16648E-02 0.25877E-02 0.58527E-02 0.53674E-02 0.65943E-02 0.74746E-02 0.80714E-02 0.82702E-02 

3 0.20250E-02 O.Z1657E-02 0.25873E-02 0.32059E-02 0.40575E-02 0.50798E-02 0.59108E-D2 0.65088E-02 0.69153E-02 0.70507E-02 

2 0.32352E-02 0.33053E-02 0.35153E-02 0.38262E-02 0.42594E-02 0.47808E-02 0.52029E-02 0.55076E-02 0.57163E-02 0.57858E-02 

I 0.44172E-02 0.44190E-02 0.44245E-02 0.44275E-02 0.44671E-02 0.45505E-02 0.45817E-02 0.45801E-02 0.45863E-02 0.45884E-02 


SOLUTION AT LZ= 4 4 0.00000 

0.00000 O.lllll 0.22222 0.33333 0.44444 0.55556 0.66667 0.77778 0.88889 1.00000 

0.83267D-16 0.38638 0.70867 1.0186 1.3661 1.7755 2.1230 2.4329 2.7552 3.1416 

10-0.68654E-03-0.18027E-03 0.13361E-0Z 0.34925E-02 0.63766E-02 0.97245E-02 0.12320E-01 0.14148E-01 0.15375E-01 0.15783E-01 

9-0.5230ZE-03-0.28365E-04 0.14525E-0Z 0.35619E-02 0.63913E-02 0.96785E-02 0.12231E-01 0.14033E-01 0.15244E-Q1 0.15646E-01 

8-0.30049E-04 0.43299E-03 0.18178E-02 0.57971E'02 0.64642E-Q2 0.95680E-02 0.11987E-01 0.13701E-01 0.14853E-01 0.15236E-01 

7 0.72588E-03 Q. 11420E-02 0.23846E-02 0.41683E-'02 0.65856E-02 Q.94047E-02 0.11612E-01 0.13183E-01 0.14241E-01 0.1<i592E-Ql 

6 0.16863E-02 0.20440E-02 0.31104E-02 0.46483E-02 0.67465E-02 0.91990E-02 0.11129E-01 0.12511E-01 0.13442E-01 0.13750E-01 

5 0.27986E-02 0.309016-02 0.39564E-02 0.52120E-02 0.69384E-02 0. 896046-02 0.10559E-(ll 0.11711E-01 0.12489E-01 0.12746E-01 

4 0.40148E-02 0.4Z349E-02 0.48865E-02 0.58353E'02 0.71528E-02 0.86977E-0Z 0.99231E-02 0.10814E-01 0.11416E-01 0.11614E-01 

3 0.52887E-02 0.54352E-02 0.58651E-02 0.64943E-02 0.73813E-02 0.84199E-02 0.92428E-02 0.98482E-0Z 0.10259E-01 0.10394E-01 

2 0.65706E-02 0.66436E-02 0.68531E-02 0.71618E-'02 0.76136E-02 0.81412E-02 0.855Z6E-02 0.88610E-02 0.90726E-02 0.91407E-02 

1 0.78237E-02 0.78358E-02 0.78674E-02 0.78968E-02 0.79476E-02 0.80192E-02 0.80163E-02 0.79691E-02 0.79364E-02 0.79241E-02 
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The output for MODE=I consists of the geometric summary information shown for 
M0DE=3. Also presented at each streamwise station are the cross sectional area of 
the duct and the mass flux. Two-dimensional tabulations of several parameters are 
also provided in the cross plane. These are: 

RADIUS - d is tance/YZERO from duct center line 
THETA - angle from symmetry plane 

VEL-S - streamwise velocity /UZERO 

VEL-TH - aximutlial velocity component/UZERO in the direction of the 
local azimuthal coordinate 
VEL-R - radial vo 1 oc i ty /UZERO 

RHO - dens i ty/RZERO 

CP/2 - half the pressure coefficient 

VOR - C(^mponent of the vorticity vector parallel to the c-enteriine 

MUT - turbulent v iscos i ty/LAMTNAR VISCOSITY 

LEN - mixing length used in turbulence viscosity calculation 

The (uitinit is presented for LEVEL 1, which is tlie results of the calculation. 
The only exception is for station 2 where the output is first presented at LEVEL 3, 
then at LEVEL 1. The LEVEL 3 output is the station I input which has been revised 
to be consistent with the governing equati(')ns. 

The logical file name convention used by the code is; 

5 Card image input 

6 Printed output 

8 Plot file NASA LeRC format 

9 Scratch storage - Z array and PFLOW 

10 Geometry for PELOW 

11 Restart 

12 Restart 

13 Pressure field 

The assignment and cataloging of these files is machine and site dependant 
and it is the responsibility of the user to make them available to the code. 
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